Characterization of topological states on a lattice with Chern number 
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We study Chern numbers to characterize the ground state of strongly interacting systems on a 
lattice. This method allows us to perform a numerical characterization of bosonic fractional quantum 
Hall (FQH) states on a lattice where conventional overlap calculation with known continuum case 
such as Laughlin state, breaks down due to the lattice structure or dipole-dipole interaction. The 
non- vanishing Chern number indicates the existence of a topological order in the degenerate ground 
state manifold. 
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One of the most dramatic manifestations of interac- 
tions in many-body systems is the appearance of new 
quantum states of matter. Often such states can be 
characterized by an order parameter and spontaneous 
breaking of some global symmetries. The 'smoking gun' 
evidence of such states is the appearance of Goldstonc 
modes of spontaneously broken symmetries. However, 
some types of quantum many-body phases can not be 
characterized by a local order parameter. Examples can 
be found in FQH systems lattice gauge theories |2j, 
and spin liquid states Q. Such states can be charac- 
terized by the topological order [J| which encompasses 
global geometrical properties such as ground state de- 
generacies on non-trivial manifolds @. Topologically or- 
dered states often exhibit fractional excitations @ and 
have been proposed as a basis of a new approach to quan- 
tum computations (7|. However, in many cases, identify- 
ing a topologically ordered state is a challenging task 
even for theoretical analysis. Given an exact wavefunc- 
tion of the ground state in a finite system, how can one 
tell whether it describes a FQH phase of a 2D electron 
gas or a spin liquid phase on a lattice? One promising 
direction to identifying topological order is based on the 
Chern number calculations[8[. The idea of this approach 
is to relate the topological order to the geometrical phase 
of the many-body wavefunction under the change of the 
boundary conditions [9( . 

Important work of Berry [Toj and Simon (Tl| initi- 
ated the investigation on geometrical phase factors and 
since then the field has been extensively studied in dif- 
ferent contexts - for a review see, for example, In 
quantum Hall (QH) systems, early works on the Chern 
number analysis [13| is focused on the Hall conductance 
and robustness of QH states against changes in the band 
structure [3] and the presence of disorder 0j|. Cur- 
rently, there is also considerable interest in understanding 
FQH states in the presence of a strong periodic potential. 
Such systems are important in several contexts including 
anyonic spin states [Ty] , vortex liquid states [17f . and ul- 
tracold atoms in optical lattices [l8l [l^. I20L l2ll. |22| which 
are promising candidates for an experimental realization. 

In this letter, we investigate a novel procedure for 
calculating Chern numbers and demonstrate that this 



method provides insight into the topological order of the 
ground state in regimes where other methods fails to pro- 
vide a definite answer for the nature of the ground state 
wavefunction. In particular, we study a fractional quan- 
tum Hall system with bosons on a lattice with a filling 
factor of v = 1/2, where v is the ratio of the number of 
magnetic flux quanta to the number of particles. In the 
continuum limit, where the flux- fraction through each 
plaquette a is very small (a <C 1), this system is ex- 
actly described by the Laughlin wavefunction. However, 
in a recent study [l8j . it was found that for some values 
of a > 0.25, the Laughlin wavefunction ceases to be a 
good description of the system, indicated by a decreas- 
ing overlap between the ground state and the Laughlin 
wavefunction. From this study, it is unclear whether this 
represented a change in the nature of the ground state, 
or just that the lattice structure distorts the state. Here, 
we use the Chern number calculation to provide an un- 
ambiguous characterization of the ground state even out- 
side the regime where there is a significant overlap with 
the Laughlin wavefunction. In particular, we show that 
the Chern number and hence the topological order of the 
system remains undisturbed until a < 0.4.(Tab|J). 

To introduce the basic ideas of the Chern number cal- 
culation in many-body system, we first discuss the con- 
tinuum regime. We consider a single particle with charge 
q on a torus T 2 (L X xL 9 ) in the presence of a magnetic 
field B perpendicular to the surface. The corresponding 
Hamiltonian is invariant under the magnetic translation 
of single particle s, i s (a) = e la ' k l h where a is a vector on 
the torus, and A; s is the pseudo-momentum of particle s, 
defined by k x{y) = -ih-g^ - qA x{y) =F qBy(x) in x and 

(y) direction, respectively, and A is the corresponding 
vector potential. The generalized boundary condition is 
given by the translation, t s (Lii)ip(x s ,y s ) = e tSi ip(x s , y s ), 
where (i = 1, 2) refer to two directions (x,y) on the torus 
T 2 and the O^s are twist angles of the boundary (FigOJi). 
The magnetic phase through each plaquette (27ra) arises 
from the field perpendicular to the surface of the torus. 
The Chern number for non-degenerate state a is defined 
by, 
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the corresponding vector field is not defined for certain 
angles (0™ , 6%) in S n regions (FigJTJj), then we should in- 
troduce a new well-defined vector field A'j (61,62), in- 
side those regions. These two vector fields differ from 



each other by a gauge transformation, A) 
\i{ a )i 







Aj (61,62) — djx(6i,6 2 ) and the Chern number re- 
duces to the winding number of the gauge transformation 
61,62) over small loops encircling (#1,6*2), i-e. dS n 
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For the case of degenerate ground state a general- 
ization of the above argument can be made, where in- 
stead of having a single vector field Aj (61,62), a ten- 



sor field A 



, 62) should be defined, a, (3 



for a g-fold degenerate ground state: A 
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Therefore, when Aj' is not defined, similar to the 
non-degenerate case, a new gauge convention should 
be acquired for those regions with singularities. This 
gives rise to a tensor gauge transformation on the bor- 



der of these regions, djX 
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h,62) and consequently the Chern number is 



given by the trace of the tensor 



FIG. 1: (a) Twist angles of the boundary condition as the re- 
sult of two magnetic fluxes threaded the torus (b) Redefining 
the vector potential around the singularities: Aj is not well- 
defined everywhere on the torus of the boundary condition. 
Therefore, another vector field A'j with different definition 
should be introduced around each singularity (9f, 9%) of Aj. 
However, Aj and A'j are related to each other with a gauge 
transformation \ and the Chern number depends only on the 
loop integrals of \ around those singularity regions. 



C{a) = — dSi d6 2 (d 1 A\ a) - d 2 A\ a> ) (1) 
27T Jo Jo 

where Aj (61,62) is defined as a vector field based 
on the eigenstate ^ a \6 l ,6 2 ) on T 2 by A < f t) (6 1 ,6 2 ) = 

In the context of QH systems, the time derivative of 
twist angles could be considered as voltage drops across 
the Hall device in two dimensions and the boundary av- 
eraged Hall conductance of the any state is related to the 
Chern number of that state er^ = C(a)e 2 /h. 

The non-trivial behavior (non-zero conductance in the 
case of quantum Hall system) occurs because of singulari- 
ties of the vector field. If for a given non-degenerate state 



C(l,2,...,q)=Y,Y I VTr X {a ^ • dt . (3) 

We focus on a system of bosons on a square lattice 
described by the Hamiltonian (l9j : 

x,y 

+ U/,n x , y (h Xt y - 1), (4) 

where J is the hopping energy between two neighbor- 
ing sites, U is the on-site interaction energy, and 2ira is 
the phase acquired by a particle going around a plaque- 
tte. We concentrate on the hardcore limit (U 3> J) and 
v = 1/2 where v is the ratio of the number of particles to 
the total number of flux in the system. The experimen- 
tal proposal for realizing such a Hamiltonian for atomic 
gases confined in an optical lattice has already been in- 
vestigated [HI, [l!| . The ground state of the system for 
very dilute lattice a < 0.2 is two-fold degenerate and is 
well described by Laughlin state. When a increases the 
lattice structure becomes more apparent and the overlap 
with Laughlin wavefunction breaks down. However, by 
numerical calculation, we show that Chern number char- 
acterizes system better and remains the same, i.e. 1/2 
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FIG. 2: Low-lying energy levels as a function of twist angles. 
For finite a the ground state energy oscillates as a function 
of twist angles and for high a > 0.4 the oscillations reach the 
exited levels, (a) shows lowest energy levels for a = 0.32 (4 
atoms on a 5x5 lattice) while (b) shows the five lowest energy 
levels for a = 0.4 ( 5 atoms on a 5x5 lattice) In both plots, 
O2 — tv and 6\ is varied from zero to 27T. 



for each state in the ground state manifold, for systems 
with higher flux density a < 0.4. 

In the case of a very dilute lattice a < 1, i.e. the 
magnetic length is much larger than the lattice spacing, 
and hence the lattice Hamiltonian approaches the con- 
tinuum limit. According to Haldane [23| . the magnetic 
translational symmetry of the center of mass results in 
a two-fold degeneracy of the ground state for v = 1/2. 
However, by increasing the magnetic field, the lattice 
structure becomes more pronounced, in such a way that 
even for a single particle, the lattice modifies the energy 
levels from being simple Landau levels into the fractal 
structure known as the Hofstadter butterfly [24| • For the 
many-body problem, the presence of the lattice causes 
the energy levels to oscillate (see Fig. [5^) and instead of 
having a unique degeneracy for all twist angles values, the 
ground state is two-fold degenerate at only certain twist 
angles. However, the two-dimensional ground state man- 
ifold is well defined and separated from the other states. 
By integrating over twist angles, one state mixes with 
the other state when levels touch each other, therefore 
in the Chern number evaluation, both levels participate 
and the degenerate form of the Chern number should be 
used (Eq. |3J). 

It is important to note that the degeneracy in the non- 
interacting regime (Landau degeneracy) is fundamentally 
different from that of the interacting hard-core case. In 
the non-interacting limit ([/ « J), the degeneracy can be 
lifted by a local perturbation e.g. an impurity, while in 
the hardcore case, the degeneracy remains in the thermo- 
dynamic limit Q . The latter degeneracy is a consequence 
of the global non-trivial properties of the manifold on 
which the particles move rather than symmetries of the 
Hamiltonian (e.g. the Ising model) [1]. Recently, it was 
shown [f| that in presence of a gap, there is a direct con- 
nection between the fractionalization and the topological 
degeneracy. In particular, the amount of the degeneracy 
is related to the statistics of the fractionalized quasipar- 



ticles e.g. in the case of v = 1/2, the two-fold degeneracy 
is related to 1 /2 anyonic statistics of the corresponding 
quasiparticles. 

The ground state degeneracy prevents the direct inte- 
gration of Eq. |T]) since wavefunctions would mix together 
when twist angles vary. Therefore, one has to use Eq. © 
and also resolve the extra gauge related to the ground 
state. We can consider two possibilities: fixing the rel- 
ative phase between the two states in the ground state, 
or lifting the degeneracy by adding some impurities. In 
the latter case, we can show that the system has a topo- 
logical order in spite of poor overlap with the Laughlin 
state [l9| . On the other hand, a significant amount of 
impurity in the system may distort the energy spectrum, 
so that the underlying physical properties of the lattice 
and fluxes could be confounded by the artifacts due to 
the impurities, especially for large a. Therefore, in this 
letter we focus on the degenerate case. 

We start with the simple case of a non-degenerate 
ground state on a discrete s-dimensional Hilbert space, 
^(01 j $2) = (ci, C2, c s ). The one-dimensional gauge 
can be resolved by making two conventions: in one 
convention the first element and in the other the sec- 
ond element of the wavefunction in the Hilbert space 
should be real i.e. we transform the ground state \& 
into = P$ = **t$ where $ = (1,0, ...,0) t is 
a s-dimensional vector and P is a projection into the 
ground state and similarly with the other reference vec- 
tor $' = (0,1, ...,0)+. Hence, we can uniquely determine 
the gauge \ which relates the two corresponding vector 
fields: e lx — <&tp$'. Therefore, the Chern number will 
be equal to the number of vorticities of x around regions 
where = <&tp<I> = |ci| 2 is zero. 

For fixing the q-dimensional ground state manifold 
gauge, we take two reference multiplets $ and 4>' which 
are two s x q matrices (q = 2 in our case). We define 
an overlap matrix as A^ = $ T P<E>, and consider the re- 
gions where det A$ or det A^, vanishes (similar to zeros 
of the wave function in the non-degenerate case). Hence, 
the Chern number for q degenerate states will be equal to 
the total winding number of Tr ^ a '^ for small neighbor- 
hoods S n , in which detA$ vanishes. It should be noted 
that the zeros of detA$ and detA^, should not coincide 
in order to uniquely determine the total vorticity. In our 
numerical calculation, we choose multiplets $ and $' to 
be two sets of two degenerate ground states at two dif- 
ferent twist angles far apart e.g. (0,0) and ("71", tt). In 
Fig. 02 we have plotted ft = det($ t P$') J detA$, and 
det A' $ , found by numerical diagonalization of the Hamil- 
tonian over a grid (30 x 30) of twist angles 6\ and 82- 
The Chern number can be determined by counting the 
number of vortices and it is readily seen that the wind- 
ing number is equal to one for the corresponding zeros of 
detA' $ and detA$. 

We have calculated the Chern number for fixed 1/ — 1/2 
and different a's by the method described above. The re- 
sult is shown in Tab HI For a C 1, we know from previ- 
ous calculation [l8[ that the ground state is the Laughlin 
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FIG. 3: (a) £1(61,62) for fixed $ and 61 and 6» 2 changes 
form zero to 2-k. This plot has been produced for 4 atoms 
in the hard-core limit on a 5x5 lattice (a = 0.32). (b) and 
(c): surface plots of detA$ and detA^ (blue is lower than 
red). 61 and 62 changes form zero to 2n. The total vorticity 
corresponding to each of trial function ($ or <&') indicates a 
Chern number equal to one for the two dimensional ground 
state manifold. 



state and we expect to obtain a Chern number equal to 
1/2 for each state i.e. total Chern number equal to one. 
For higher a, the lattice structure becomes more appar- 
ent and the overlap with the Laughlin state decreases 
rapidly. However, in our calculation, the ground state 
remains two-fold degenerate and the associated Chern 
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TABLE I: Chern Number for different configurations in the 
hard-core limit for fixed filling factor v = 1/2. The Laughlin 
state overlap is shown in the last column, although it deviates 
from the Laughlin state. Although the ground state deviates 
from the Laughlin state, the Chern number remains equal to 
one half per state before reaching some critical a c ~ 0.4 where 
the energy gap vanishes. 



number remains equal to one before reaching some criti- 
cal a c ~ 0.4. Hence, we expect to have similar topolog- 
ical order and fractional statistics of the excitations on 
the lattice in this regime. 

For higher flux densities, a > a c , the two-fold ground 
state degeneracy is no longer valid everywhere on the 
torus of the boundary condition. In this regime, the is- 
sue of degeneracy is more subtle, and finite size effect 
becomes significant. The translational symmetry argu- 
ment [23| is no longer valid and the degeneracy of the 
ground state varies periodically with the system size [25[ . 
Some gaps might be due to the finite size and vanish in 
the thermodynamic limit. To investigate this, we study 
the ground state degeneracy as a function of boundary 
angles (61,62) which are not physical observable. There- 
fore, the degeneracy in thermodynamic limit should not 
depend on the their value. In particular, Fig. ^jp shows 
the energy levels of five particles at a — 0.4 for different 
values of twist angles. The first and the second level are 
split at (61 = 82 = 0), while they touch each other at 
(61 = 62 = 7t). We have observed similar behavior for 
different number of particles and lattice sizes e.g. 3 and 
4 atoms at a = 0.5. Therefore, the ground state enters a 
different regime which is a subject for further investiga- 
tion. Existence of the topological order does not require 
a very strong interaction i.e. hard-core limit. Even at 
finite interaction strength U ~ J a, we have observed the 
same topological order with the help of the Chern num- 
ber calculation. If U gets further smaller, the energy gap 
above the ground state diminishes and the topologi- 
cal order disappears. 

One of the impediment of the experimental realization 
of Quantum Hall state is the smallness of the gap which 
can be improved in presence of the dipole-dipole inter- 
action (l9| . The dipolc interaction can be represented 
as extra term . Ud ru nj/|i"i — rj| 3 in the Hamiltonian 
Eq. ([4]), where is the number of particles at location 
ri in the units of lattice spacing and Ud is the strength 
of the interaction. The magnetic dipole-dipole interac- 
tion has been achieved in Bose-Einstein condensation of 
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Chromium 26], however, for a lattice realization, polar 
molecules with strong permanent electric dipole moments 
are more promising candidates, where the dipole interac- 
tion can be an order of magnitude greater than the tun- 
neling energy. In the presence of such strong long-range 
interaction, the ground state deviates from the conven- 
tional FQH state even in the continuum case (i.e. for 
even a < 0.2, the overlap with the Laughlin wavefunc- 
tion decreases by increasing the strength of the dipole 
interaction). However, by evaluating Chern number, we 
are able to identify the topological order of the system, 
that turns out to be intact i.e. Chern number equal to 
one for the two-fold degenerate ground state. 

In conclusion, we have investigated a method to un- 
ambiguously calculate the Chern number for the ground 



state of a system. For the FQHE system on a lattice that 
we have investigated, the Laughlin wavefunction ceases 
to be a good description of the ground state for high 
fluxes a > 0.25, but the Chern number remains 1/2 per 
state until a < 0.4 which is a direct indication of topo- 
logical order in the system. Calculating Chern numbers 
by this method can be generalized for finite lattice sys- 
tems to properly characterize the ground state manifold 
which is otherwise impossible with conventional overlap 
methods. 
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